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Abstract. We calculate a grid of star models with and without the effects of axial rotation for stars in the mass 
range between 2 and 60 Mq for the metallicity Z = 10~^. Star models with initial masses superior or equal to 
9 M0 were computed up to the end of the carbon-burning phase. Star models with masses between 2 and 7 
M0 were evolved beyond the end of the He-burning phase through a few thermal pulses during the AGB phase. 
Compared to models ai Z = 0.02, the low Z models show faster rotating cores and stronger internal fi-gradients, 
which favour an important mixing of the chemical elements. The enhancement of N/C at the surface may reach 2 
to 3 orders of magnitude for fast rotating stars. Surface enrichments may make the evolved stars less metal poor 
than they were initially. In very low Z models, primary nitrogen is produced during the He-burning phase by 
rotational diffusion of ^^C into the H-burning shell. A large fraction of the primary ^'*N escapes further destruction 
and enters the envelope of AGB stars, being ejected during the TP-AGB phase and the formation of a planetary 
nebula. The intermediate mass stars of very low Z are the main producers of primary ^*N, but massive stars also 
contribute to this production; no significant primary nitrogen is made in models at metallicity Z=0.004 or above. 
We calculate the chemical yields in He, C, N, O and heavy elements and discuss the chemical evolution of the 
CNO elements at very low Z. Remarkably, the C/O vs 0/H diagram is mainly sensitive to the interval of stellar 
masses, while the N/O vs O/H diagram is mainly sensitive to the average rotation of the stars contributing to 
the element synthesis. The presently available observations in these diagrams seem to favour contributions either 
from stars down to about 2 Mq with normal rotation velocities or from stars above 8 M0 but with very fast 
rotation. 
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1. Introduction 

Stellar rotation really modifies all the o utputs of stellar 
evoluti on for massive stars (Hegeret al. 2000a ; Heger & 
Langer |2000l Meynet & Maeder ^OOOl ). This is true at so- 
lar metallicity. At lower metallicitics, like Z=0.004 in the 
SMC, we noticed that the effects of ro tation are expected 
to be larger (Maeder & Meynet 2001 ). In particular, for 
similar initial distributions of the rotational velocities, a 
larger fraction of the stars at lower Z reach break-up ve- 
locities. This is a result of the smaller losses of angular 
momentum by the stellar winds. 

In addition, it may be that the initial distribution of 
the rotational velocities is not the same at lower Z. Indeed, 



it has been shown (Maeder et al. [1999D that the fraction of 
Be stars (i.e. stars close to break-up) is much higher in the 
SMC than in the Milky Way. However, we do not know 
whether this is just a consequence of the smaller mass 



loss, as said above, or whether the initial distribution of 
rotation velocities is also different as a result of processes 
of star formation at low Z. Whatever the exact origin of 
the higher fraction of stars close to break-up at lower Z, 
this shows the need of studies of star models with rotation 
at low Z. 

We consider here star models with metallicity Z — 
10~^, which is low enough to correspond to the most ex- 
treme metallicity observed in halo stars of the order of 
[Fe/II]~ -3 and which nevertheless avoid the particular- 
ities oi Z — models, which we may consider in a fu- 
ture paper. The main possible comparisons with the ob- 
servations will concern the chemical evolution of the abun- 
dances of the CNO elements and of other heavy elements 
in halo stars and very low Z galaxies. This is why we 
put here a particular emphasis on the chemical yields in 
CNO at very low [Fe/H]. This is a topical point in rela- 
tion with the problem of primary nitrogen (Edmunds and 
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Fig. 1. Stream lines of meridional circulation in a rotating 20 Mq model with solar metallicity and v-mi = 300 km 
at the beginning of the H-burning phase (see text). The streamlines are in the meridian plane. In the upper hemisphere 
on the right section, matter is turning counterclockwise along the outer stream line and clockwise along the inner one. 
The outer sphere is the star surface and has a radius equal to 5.2 Rq. The inner sphere is the outer boundary of the 
convective core. It has a radius of 1.7 R©. 



Pa gel |1978| ; Barbuy |1983| ; Carb on et al. |l987t Th uan et 



al 



1995 



Izotov & Thuan |1999[ Henry et al.|00^. Also, 
the recent debate around the behavior of the [0/Fe] ra- 
tio at very low Z (cf. Israelian et al. 2001; Melendez et 
al. 2001) shows the need of a better understanding of the 
CNO yields at very low metallicities. 

In Sect. 2, we discuss the model physics. In Sect. 3, 
we examine the internal rotation and the surface veloci- 
ties in Sect. 4. The models with zero rotation are briefly 



mentioned in Sect. 5. The HR diagram and lifetimes are 
discussed in Sect. 6. The evolution of surface abundances 
are examined in Sect. 7. In Sect. 8, we discuss the prob- 
lem of the origin of primary nitrogen and we show how 
rotation can solve it. The chemical yields in He, CNO and 
heavy elements are discussed in Sect. 9. 
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Table 1. Initial abundances in mass fraction. 



Element 


Initial abundance 






3 He 




•*He 


0.22999812 


12(~, 


7.5500C-7 


13q 


O.lOOOc-7 




2.3358C-7 




0.0092C-7 


16o 


67.100C-7 


I'O 


0.0300C-7 




0.1500C-7 




7.8368C-7 




0.6306C-7 


'"Mg 


3.2474C-7 


^=Mg 


0.4268C-7 




0.4897C-7 



2. Physics of the models 



The initial 
composition 
paper VII, 
Rogers 
tures 



composition is 
is enhanced 
the opacities 



given in Table The 
in a-elements. As in 
are from Iglesias & 



(1996), complemented at low tempera- 
with the molecular opacities of Alexander 
( http: / /web. physics. twsu.edu/alex/wwwdra.htm| ) . The 
nuclear reaction rates are also the same as in paper VII 
an d are based on the new NACRE data basis ( Angulo et 
al. 119991) . 

The physics of the present models at Z — 10^^ is the 



same as for models at Z=0.004 (Maeder & Meynet |200l| ). 
For rotation, the hydrostatic effects and the surface dis- 
tortion are included (Meynet & Maeder 1997), so that 
the Toff given here corresponds to an average orientation 
angle. The diffusion by shears, which is the main effect 
for the mixing of chemical elements, is included (Maeder 



1997), with account of the effects of the horizontal turbu- 
lence, which reduces the shear effects in regions of steep 
/z-gradients and rei nforce s it in regions of low /x-gradients 
(Maeder & Meynet ^00l|) . 

Meridional circulation is the main effect for the in- 
ternal transport of angular momentum. We use here the 
expression by Maeder & Zahn (1998) for the vertical com- 
ponent U{r) of the meridional circulation. It is interesting 
to represent graphically this circulation. Fig. [l] illustrates 
the patterns of the meridional circulation in a 20 Mq star 
at Z = 0.020 and initial rotation velocity Vini — 300 km 
s~^ on the ZAMS. The figure is symmetrical with respect 
to the rotation axis, as well with respect to the equatorial 
plane. The small inner sphere is the edge of the central 
convective core. The inner tube, in the upper hemisphere, 
represents an interior cell of the meridional circulation. 
There is an ensemble of such concentric tubes with differ- 
ent meridional velocities. The motions occur in a meridian 
plane (i.e. turning around the tube). In the upper hemi- 
sphere and around the inner tube, the fluid elements go 
upward on the inner side of the tube and descend toward 



o 

C3 




Fig. 2. Evolution of the angular velocity as a function 
of the distance to the center in a 15 Mq star with Vi^i 
= 300 km s~^ and Z = 10~^. Xc is the hydrogen mass 
fraction at the center. The dotted line shows the profile 
when the He-core contracts at the end of the H-burning 
phase. 



the equator on the outer side of the tube (U{r) is pos- 
itive). The external tube represents an outer circulation 
cell, due to the Gratton-Opik term, which is important in 
the outer stellar layers. This term leads to a negative U{r), 
which means that, in the upper hemisphere, the fluid goes 
up on the outer side of the tube and down along the inner 
side. There also, this tube is one among an ensemble of 
stream lines turning in the meridian plane. 

The mass loss rates are based on the same references as 
in the paper for the Z=0.004 models (Maeder & Meynet 
^001 ) and in particular on the data by Kudritzki & Puis 
( ^OOOD for the OB stars. Of course, the strong reduction of 
the mass loss rates with metallicity for stars below 60 Mq 
makes the mass loss rather unimportant for the metallicity 
Z = 1Q~^ considered here, as illustrated by Table 1 which 
shows the values of the final masses. We account for the 
effects of rotation on the mass loss rates, according to the 
standard stellar wind theory applied to a rotating star 



(Maeder & Meynet gOOOj). The net result is that the very 
massive stars with initial M > 60 Mq may still experience 
significant mass loss as shown in Table 1, if they rotate 
very fast. 

Star models with initial masses superior or equal to 9 
Mq were computed up to the end of the carbon-burning 
phase. Star models with masses between 2 and 7 Mq were 
evolved beyond the end of the He-burning phase through 
a few thermal pulses during the AGB phase. 
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Fig. 3. Variation of the angular velocity as a function 
of the distance to the center in 3 and 9 star models 
with Wini = 300 km s'^ at Z = 0.020 and Z = IQ-^. The 
mass fraction of hydrogen at the centre ~ 0.40. 

3. The evolution of the internal rotation and 
meridional circulation 

There are remarkable differences in the internal distribu- 
tions of the angular velocity ri(r) depending on the stellar 
metallicity Z . This was already suggested in paper VII 
(Maeder and Meynet 2001 ), when comparing models at 
Z = 0.004 and Z = 0.020. It is extended here with models 
at Z = 10-^ 

These matters are not academic problems ! Indeed, the 
distribution of determines for example the mixing 

of chemical elements, the size of the convective core and 
therefore the chemical yields. The results in Sect. 8 and 
9 below on the chemical yields are a consequence of the 
distribution of ri(r). 

Fig. ^ shows the evolution of ^ir) during the MS phase 
of a 15 M0 a,t Z = 10"^, (this follows the initial conver- 
gence of the -profile which is very short, i.e. < 1% of the 
MS lifetime). We notice that the rotation of the convective 
core only has a small decrease during the MS phase, much 
smaller than at higher metallicities. This results from 2 
effects, a) The mass loss a,t Z = 10~^ is much smaller 
than at solar composition and thus less angular momen- 
tum is removed from the star, b) As we shall see below, 
the meridional circulation is very slow in the outer regions 
of the models at very low Z and it transports much less 
angular momentum outwards than in models at solar com- 
position. In view of these remarks, it is likely that massive 
stars at lower Z have faster spinning cores. 

Another significant difference shown by Fig. ^ concerns 
the gradient of fl outside the convective cores. Here, the 
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Fig. 4. Evolution of U{r) the radial term of the vertical 
component of the velocity of meridional circulation for a 
model of 20 Mq with Z = 10~^ at various stages during 
the MS phase. Xc is the hydrogen mass fraction at the 
center. The dashed line shows the values of U{r) inside a 
20 Mq model at Z = 0.004 when = 0.28. 



gradients are steeper and they remain significant up to 
the stellar surface, while a.t Z — 0.02 the fi-distribution 
becomes very flat in the ex terna l layers, as evolution pro- 
ceeds (Meynet & Maeder 2000 ). This difference is well 
illustrated in Fig. ^, where we notice for the 3 and 9 Mq 
models the much steeper fi-gradients at lower Z, while 
the models at Z — 0.02 show very flat gradients in the 
outer layers. The reason for the higher f2-gradients here 
are the same as for the faster spinning cores. These higher 
f2-gradients imply stronger shears and thus more mixing 
by shear diffusion, which is the main effect for the out- 
ward transport of the chemical species. (The differences 
in between the 9 and 3 Mq models result from the fact 
that we consider stars with the same Wini, but different 
radii). 

Fig. ^ shows an example at Z = 10~^ of the evolution 
of U{r), the vertical component of the velocity of merid- 
ional circulation. The size and evolution of U{r) is very 
different from the case at Z = 0.02. At Z = 0.02, U{r) 
takes large negative values particularly in the outer lay- 
ers. This is due to their low density, which makes a large 
Gratton-Opik term in the expression of U{r), (cf. 

Maeder and Zahn 199^ ). At Z = 10~^, the large negative 



values of C/(r) have disappeared, U{r) is equal to lO""^ cm 
at the end of the MS phase, while it was 50 times more 
negative in the c orresponding models at Z — 0.02 (Meynet 
& Maeder 2000 ). The differences do not concern so much 
the deep interior, but mainly the outer layers. The physi- 
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Fig. 5. Evolution of the surface equatorial velocity as a 
function of time for stars of different initial masses with 
Vini = 300 km s~^ and Z = 10~^. The track without label 
corresponds to a 20 Mq model. 

cal reason of the above differences is the fact that the star 
is more compact at lower Z and that the density in the 
outer layers is not as low as at solar composition. 

Fig. ^ also shows the curious curve for a model at Z = 
0.004. In the interior, U{r) is about the same as in the 
present models (and this is true for all Z values) . The big 
external dip of U{r), which was present at Z = 0.02 is 
very much reduced, but still present at Z — 0.004, while 
at Z — 10~^ the external dip is fully absent. 

The small U{r) in the external layers of the present 
models is mainly responsible for the presence of an f2- 
gradient up to the stellar surface (cf. Fig. Since the 
mixing of the chemical elements is mainly driven by the 
shear, the presence of this il-gradient in the outer layers 
enables the large mixing and surface chemical enrichments 
that are present in the Z = 10~^ models. 

4. The evolution of the surface rotation velocities 

It will probably be a certain time until we are able 
to observe wsini for stars in galaxies with Z = 10^^. 
Nevertheless, these objects have contributed to shape the 
composition of our universe and they deserve some inter- 
est. 

Fig. 1^ and Fig. ^ show the evolution of the surface ro- 
tational velocities for models with initial masses from 2 to 
60 Mq. Figs. I and H show the corresponding evolution of 
jy- . We notice the relative constancy of jy- during the MS 
evolution for stars with mass between 5 and 20 Mq. The 
cases of 40 and 60 Mq are noticeable as shown by Figs. || 
and 0. These models reach the break-up velocities near 
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Fig. 6. Same as Fig. || for a 2, 3 and 5 Mq stellar model. 




years] 



Fig. 7. Evolution of the ratio il/ilc of the angular velocity 
to the break-up angular velocity at the stellar surface for 
stars of different masses at Z ~ 10^^. 



the end of the MS phase. This is completely different from 
the models at Z = 0.02, where the rotation becomes very 
small due to the huge losses of mass and angular momen- 
tum. Thus, if the initial mass function at low Z extends up 
to high mass stars, as often supposed, rotation is likely to 
be a major effect in the course of the evolution of massive 
stars, since many of them are likely to reach break-up ve- 
locities. This would even more be the case for the massive 
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Fig. 8. Same as Fig. for a 2, 3 and 5 M0 stellar model. 
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At low Z like in the models at Z = 10~^ , the growth of 
is possible because of the very small mass loss and 
also it is favoured by the outward transport of angular 
momentum which is much larger for the more massive 
stars. As shown by Maeder and Meynet (2001), the values 
of U(r) are more negative for the larger stellar masses, due 
to several facts: lower gravity, higher radiation pressure, 
larger L/M ratio and especially the lower density. Thus, 
the outward transport is more efficient. 

In view of these results, we may somehow precise our 



suggestion (Maeder and Meynet 2001) that at very low Z 
a large fraction of the massive stars reach their break-up 
velocities. This seems true for the highest masses above 
about 30 M0, but not necessarily for the OB stars below 
this limit. This question is of high importance, because 
if the massive stars reach their break-up velocity, most 
of their evolutionary and structural properties will be af- 
fected. For example, they could also lose a lot of mass 
and produce some Wolf-Rayet (WR) stars. They would 
have a relatively small remaining mass at the time of the 
supernova explosion, like their counterparts at solar com- 
position. 

For the models of 3 Mq, as illustrated in Figs. || and 
^, the rotation velocity remains about constant during the 
MS phase, while ^ increases. For 2 Mq, we notice a net 
increase. This particular behavior is due to the different 
shape of the track of the 2 Mq model in the HR diagram 
(Fig. |o|) , which mimics the tracks of lower masses dom- 
inated by the pp chain. This is well explainable, because 
at Z = 10^'"' the CNO cycle is less important than at so- 
lar composition, thus the mass limit where the CNO cycle 
starts dominating over the pp chain is shifted upward. For 
this model of 2 Mq , there is no large increase of the radius 
during the MS evolution and thus rotation keeps higher. 

When one examines the evolution of v sin i during the 
MS phase for stars of the same mass but different initial 
velocities, one usually notes at Z— 0.02 a so-called velocity 
convergence (Langer 1998| ). This is due to the fact that 
the faster rotating stars lose more mass and thus more 
angular momentum. In the present models, the mass loss 
rates are in general very small (as long as the stars are not 
at break-up) and thus there is no velocity convergence, i.e. 
the stars of different initial velocities finish the MS phase 
with different velocities as illustrated by Table 1. 



Fig. 9. Evolution of the surface equatorial velocity as a 



function of time for 60 Mq stars with Vini 
at different initial metallicities. 



300 km s" 



stars which have a blueward evolution as a result of strong 
internal mixing. Their radii would decrease, thus favour- 
ing extreme rotation velocities. We note that for masses 
between 3 and 20 Mq, the rotation velocity keeps about 
constant during the MS phase, before decreasing in the 
post-MS phases. 

Fig. ^ clearly illustrates the very different evolution of 
the rotational velocities of a 60 Mq at various metallicities. 



5. Models with zero rotation 

For purpose of comparison, we have computed non- 
rotating stellar models with the same physical ingredients 
as for the computation of the rotating ones. The evolu- 
tionary tracks are presented in Fig. |l^, the lifetimes and 
various properties of the models are given in Table § The 
models were computed with the Schwarzschild criterion 
for convection and therefore present the usual differences 
when compared with models accounting for the effect of 
overshooting (see paper V and Vll for a more detailed 
discussion). 
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Fig. 10. Evolutionary tracks for non-rotating (dotted lines) and rotating (continuous lines) models for a metallicity 
Z = 10^^. The rotating models have an initial velocity Vini of 300 km s^^. 



With respect to our previous grids of stellar models 
at solar metallicity (see paper V), the present models 
are shifted toward higher effective temperatures, by about 
0.15-0.20 dex in log T^g. The stars are much more com- 
pact than at solar metallicity by nearly a factor two (more 
precisely by a factor between 1.8-1.9 depending on the 
initial mass). This is a well known consequence of the low 
opacities in the outer layers of metal poor stars. 

In paper VII, we noticed that models at Z — 0.004 
with initial masses between 10 and 12.5 were showing 
a behavior in the HR diagram intermediate between the 
cases of stars presenting a well developed blue loop and the 
case of more massive stars, which do not produce any blue 
loop, but begin to burn their helium in their core at a high 
effective temperature, while they cross the HR diagram 
for the first time. Here aX Z = 10~^, this transition zone 
covers a broader range of initial masses, from 2 to about 
15 Mq. This is well consistent with grids of stellar models 
computed by other authors. Indeed a similar broadening 
in mass for this transition region can be observed in the 
grids by the Padova group (compare for instance the grids 



by Fa gotto et al. |1994| at Z = 0.0004 and that of Girardi 
et al. |199^ at Z = 0.0001). 

From Fig. |o|, one sees that the 20 Mq model does not 
reach the red supergiant phase at least before the end of 
the C-burning phase. This feature is also present in zero 
metallicity stellar models (see e.g. Marigo et al. 2001), but 
for pop III models, it extends over a broader range of ini- 
tial masses (from ~ 10 to ~ 40 Mq). Finally, let us note 
that, when the metallicity decreases, the mass limit for 
the helium flash decreases (see also Marigo et al. 2001). 



This is a consequence of the higher central temperatures 
reached in metal poor stars. Typically a 2 Mq model at 
solar metallicity, computed with the same physical ingre- 
dients as used in the present work, undergoes He-flash, 
while the corresponding model at Z = 10^^ ignites he- 
lium in a non-degenerate environment. 

6. HR diagram, mass-luminosity relations and 
lifetimes 



The effects of rotation a t Z — 0.020 have already be en dis 
cussed in Talon et al. ( |l997|) , Denissenkov et al. ( 1999 ) 
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Table 2. Properties of the stellar models at the end of the H-burning phase, at the end of the He-burning phase 
and at the end of the C-burning phase or during the thermal pulse-AGB phase. The masses are in solar mass, the 
velocities, in km s^^ and the lifetimes, in million years. The abundances are in mass fraction. The abundance ratios 
are normalized to their initial values, which are, in mass fraction, (N/C)ini== 0.309 and (N/0)ini= 0.035 . 
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Models at the beginning of the C— burning phase. 



Heger et al. (2000a| ), Heger & Langer ( |2000| ), Meynet & 
Maeder ( |200C ). At the metallicity of the SmaU Magellanic 
Cloud, the effects of rotation have been discussed by 
Maeder & Meynet ( ^00 1] ). Let us very briefly recall the 
most important effects of rotation: 

— The Main Sequence (MS) evolutionary tracks with ro- 
tation are extended toward lower effective tempera- 
tures and reach higher luminosities at the end of the 
MS phase than their non-rotating counterparts. In 
that respect rotation acts as a moderate overshoot. 

— The MS lifetimes increase with rotation, typically by 
about 10% for a 2OM0 model with an average rota- 
tional velocity on the MS corresponding to the ob- 
served ones. 

— For a given value of the initial mass and metallicity, 
the evolutionary tracks may be different, due to differ- 
ent initial rotational velocities. Even for a given initial 
rotation, the tracks may appear different depending on 
the angle of view, since the polar regions of a rotating 
star are in general hott er tha n the equatorial ones (cf. 
Maeder & Peytremann 197C ) . This induces some scat- 



ter in the position of the end of the MS phase in the 
HR diagram. 

— The theoretical period-luminosity relation for the 
Cepheids is changed by rotation. A Ccpheid, at a given 
position in the HR diagram, if originating from a rotat- 
ing progenitor, will pulsate with a longer period than 
a Cepheid having a non-rotating progenitor. 

— The evolution toward the red supergiant phase is 
favoured by rotation. Rotation also facilitates the for- 
mation of Wolf-Rayet stars. 

— The surface abundances are modified. 



We shall see that the models at the very low metallicity 
Z = 10^^ present some striking differences with respect to 
what happens at higher metallicities: in particular, rota- 
tion implies smaller main sequence lifetimes and rotation 
does not favour the evolution toward the red supergiant 
stage (at least for the range of initial velocities explored 
here). 
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6.1. The HR diagram 

For most of the stellar models, we computed the rotating 
tracks for an initial velocity Wini = 300 km . This value 
of Uini corresponds to a mean velocity v during the MS 
between 226 and 240 km s~'^ for initial masses below 20 
M0 (see Table ||). These values are close to the mean 
rotational velocities observed for OBV type stars at solar 
metallicity, which are between 200-250 km s~^. For more 
massive stellar models aX Z — 10~^ the average velocities 
are higher. This results essentially from the larger outward 
transport of angular momentum by circulation i n mor e 
massive stars (see Section 4 and Maeder & Meynet 2001 ). 

Figure |l^ shows the evolutionary tracks of non- 
rotating and rotating stellar models for initial masses be- 
tween 2 and 60 Mq. The effective temperatures plotted 
correspond to an average orientation angle (see also pa- 
per I). At the beginning of the evolution on the ZAMS, 
rotational mixing has no impact on the structure, since 
the star is homogeneous. At this stage, only the hydro- 
static effects of rotation are present, i.e. the effects due to 
the centrifugal acceleration term in the hydrostatic equi- 
librium equation. As is well known, these effects shift the 
ZAMS position toward lower values of L and Toff (see e.g. 
paper I). From Fig. one sees that the less massive the 
star, the greater the shift. This results from the facts that, 
for a given Wini, the lower the initial mass, the greater the 
ratio of the centrifugal force to the gravity. Indeed this ra- 
tio, equal to ^^vJ, varies as about 1/Af" with a equal 
to about 0.4. 

As was the case at higher metallicities, the MS width 
is increased by rotation. Rotational mixing brings fresh 
H-fuel into the convective core, slowing down its decrease 
in mass during the MS. A more massive He-core is pro- 
duced at the end of the H-burning phase, which favours 
the extension of the tracks toward lower effective tem- 
peratures. Rotational mixing also transports helium and 
other H-burning products (essentially nitrogen) into the 
radiative envelope. The He-enrichment lowers the opacity. 
This contributes to the more rapid increase of the stellar 
luminosity during the MS phase and limits the redwards 
motion in the HR diagram. 

The widening of the MS produced by rotation mimics 
the effect of an ov ershoot beyond the convective core (see 
Talon et al. 1997 , paper VII). Since the observed width of 
the MS has often been taken to parameterize the size of the 
convective core, the above argument shows that rotating 
models tend to decrease the amplitude of the overshoot 
necessary to reproduce the observed MS width. 
Figure 
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shows the evolutionary tracks of 20 Mq mod- 
els for different initial velocities and metallicities during 
the H-burning phase. One sees that the extension of the 
MS due to rotation decreases when the metallicity de- 
creases. This results from the smaller increase of the He- 
core due to rotation at low Z. Typically, at the end of 
the MS at Z = 10~^, the helium core mass in the rotat- 
ing 20 Mq model {vini = 300 km s^^) is greater by 11% 
with respect to its value in the non-rotating model. The 



corresponding increase at Z=0.004 is 23%. The reason for 
this difference is the following one. When rotation brings 
fresh hydrogen fuel, in the core, it brings also carbon and 
oxygen which act as catalysts in the CNO burning. These 
catalyst elements are of course in much lower abundances 
dX Z — 10^^ than at Z = 0.004 and thus the core en- 
hancement is less pronounced. 

At Z — 0.004, the rotating star models with initial 
masses between 9 and ~25 M0, are evolving, after the 
MS phase, much more rapidly toward the red supergiant 
stage (RSG) than the non-rotating models, in agreement 
with the observed number ratio of blue to red supergiants 
in the Sma ll Ma gellanic Cloud cluster NGC 330 (Maeder 
& Meynet 2001). As was extensively discussed, the bal- 



ance between the blue and the red is very sensitive to 
many effects. At the very low metallicity considered here, 
rotation does not succeed in producing red supergiants, at 
least for the range of initial rotational velocities explored. 
The dominant reason appears to be the smaller He-cores 
produced at lower Z and the fact that the growth of this 
core due to rotation is smaller than at higher Z . In terms 
of the discussion by Maeder & Meynet ( ^001 ) , this reduces 
the central potential enough to keep a blue location during 
the whole He-burning phase, and the amount of helium 
diffused in the region of the shell is unable to compen- 
sate for the smaller core. Only, during the very last stages 
in the C-burning phase, when the core heavily contracts, 
does the central potential grow enough to produce a red 
supergiant. 

The rotating models do not present any well devel- 
oped blue loop, except in the case of the 2 Mq model. 
In that respect the situation is similar to the case of the 
non-rotating models (see Sect. 5 above). Interestingly, we 
note that the rotating models, with initial mass below 7 
M0, evolve redwards during the AGB-phase. This a conse- 
quence of the third dredge-up, which brings at the surface 
carbon and oxygen synthesized in the He-burning shell, as 
well as primary nitrogen built up in the H-burning shell 
(see Sect. 7 and 8). The important enhancements of these 
elements at the surface make the star to behave as a more 
metal rich star and thus push it to a redder location in 
the HR diagram. 

In the present grid no model enters the Wolf-Rayet 
phase. At the end of the C-burning phase, the mass frac- 
tion of hydrogen at the surface of the rotating 60 Mq 
model is still important ('^ 0.48), although much lower 
than at the surface of the corresponding non-rotating 
model (^ 0.77). It is likely that more massive or faster 
rotating star models may enter the Wolf-Rayet phase be- 
fore central He-exhaustion. 



6.2. Masses and mass-luminosity relations 

Table || presents some properties of the models. Columns 
1 and 2 give the initial mass and the initial velocity Vi^i re- 
spectively. The mean equatorial rotational velocity v dur- 
ing the MS phase is indicated in column 3. The H-burning 
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Fig. 11. Evolutionary tracks for rotating 20 M0 models 
with different initial velocities and various initial metal- 
licities. The initial velocities Uini are indicated. See Table 
1 for more details on the models at Z = 10^^. 



lifetimes in, the equatorial velocities u, the helium surface 
abundance Yg and the surface ratios N/C and N/0 at the 
end of the H-burning phase and normalized to their ini- 
tial values are given in columns 4 to 8. The columns 9 
to 13 present some properties of the models at the end of 
the core He-burning phase; ^hc is the He-burning lifetime. 
Some characteristics of the last computed models are given 
in columns 14 to 18; Mfjn is the final stellar mass. For stars 
with initial mass superior or equal to 9 Mq , the final stage 
corresponds to the end of the C-burning phase. For the 
lower initial mass stars, it corresponds to the beginning of 
the Thermal Pulse AGB (TP-AGB) phase. Typically the 
rotating 3 M0 model was computed until the fifth thermal 
pulse. For the intermediate mass stars, the mass fractions 
of carbon (C), oxygen (O) and nitrogen (N) at the surface 
of the stars are given. 

Rotation, by enhancing the luminosity and lowering 
the effective gravity, increases the mass loss rates (Maeder 
& Meynet 2000). As a consequence, the final masses of the 
rotating models are smaller. At the metallicity Z = 10~^, 
except for the 60 Mq model, the effects of rotation on the 
final stellar masses are very weak (see Table ||). 

In general, rotation makes the star overluminous for 
their actual masses. Typically for Vi^i = 300 km s'-^, the 
luminosity vs. mass (L/M) ratios at the end of the MS are 
increased by 10-14% for stars in the mass range from 3 
to 40 Mq. This results essentially from the He diffusion 
in the radiative envelope which lowers the opacity and 
makes the star overluminous. In the 60 Mq model, mix- 
ing is particularly efficient and the increase of the L/M 
ratio amounts to 23%. For the 2 M© model, the L/M ra- 



tio is decreased in the rotating model, by 6-7% . In this 
last case, the convective core during the H-burning phase 
disappears very early, when the mass fraction of hydrogen 
in the center is still high {X^ — 0.52 in the Vini — 300 km 
model). This puts farther away from the surface the 
region where helium is produced and thus slows down the 
helium diffusion in the outer envelope. 

The increase of L/M, due to rotation, a,t Z — 10^^ are 
in general inferior to those obtained at Z = 0.004, which 
are between 15-22% (Maeder & Meynet 2001). This is 
mainly a consequence of the following fact: at Z = 10~^, 
the increase of the H-burning convective core due to rota- 
tion is inferior to that obtained at Z ~ 0.004 for the same 
value of Willi. 

6.3. Lifetimes 

Generally we can say that the MS lifetime duration is 
affected by rotation at least through three effects: 

— 1) Rotation increases the quantity of hydrogen burnt 
in the core. This increases the MS lifetime. 

— 2) The hydrostatic effects of rotation make a star of a 
given initial mass to behave as a non-rotating star of 
a smaller initial mass. This tends to increase the MS 
lifetime. 

— 3) Rotation increases the helium abundance in the 
outer radiative envelope. This tends to make the star 
overluminous with respect to its non-rotating counter- 
part and thus to reduce the MS lifetime. 

When the metallicity decreases, the effect number 3 tends 
to become the most important one. 

This can be seen from a detailed comparison of the 
tracks in Fig. |l^. Indeed the evolutionary tracks for our 
rotating 20 M© models (wini = 300 km s~^) become overlu- 
minous with respect to the non-rotating tracks at an ear- 
lier stage for lower metallicities. This mainly results from 
the fact that when the metallicity decreases, rotational 
mixing is more efficient (see Sect. 3). Also, at lower Z, the 
stars are more compact and therefore the timescale for 
mixing, which is proportional to the square of the radius, 
decreases. This favours the helium diffusion in the outer 
envelope. The diffusion of hydrogen into the core, which 
would increase the MS lifetime, is not really affected, be- 
cause hydrogen just needs to migrate over the convective 
core boundary to be engulfed into the core. One notes 
also that for given values of the equatorial velocity and 
of the initial mass, the ratio fi/ilc of the angular velocity 
to the break-up velocity decreases whith the metallicity. 
Thus the hydrostatic effects, which usually make the star 
fainter, are in general smaller at lower Z. 

As a consequence of the above effects, at Z = 10^^, the 
MS lifetimes are decreased by about 4-14% for the mass 
range between 3 and 60 M© when Vini increases from to 
300 km s-i (cf. Table I). 

We notice that the rotating 2 M© model has a longer 
MS phase than its non-rotating counterparts, in contrast 
with what happens for higher initial mass stars. This is 
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Fig. 12. Evolution as a function of logTcff of the abun- 
dance ratio N/C where N and C are the surface abun- 
dances of nitrogen and carbon respectively. The abun- 
dance ratios are normalized to their initial values. The 
tracks are for 9 M© for different values of the metallic- 
ity, Z, and rotation. The long-dashed line, at the bot- 
tom, corresponds to a non-rotating 9 Mq stellar model at 
Z = 10"^ 

because, when the initial mass decreases, the hydrostatic 
effects become more and more important (see Fig. |lo|) . 

For what concerns the effects of rotation on the He- 
burning lifetime, let us simply say that when vini increases 
from to 300 km s~^, the changes in the He-burning 
lifetimes are inferior to 10%. 



7. Evolution of the chemical abundances at the 
surface 

Fig. |l^ shows the evolution of the N/C ratios in models 
of rotating stars with 9 M© for the initial Z = 0.02, 0.004 
and 10~^. Values for other stellar masses may be found in 
Table || 

At zero rotation, for any Z and any masses there is 
no enrichment during the MS phase (except at Z—0.02 
for M > 60Mq due to very high mass loss). At 9 M© for 
an initial rotation of 300 km s~^, we notice an increase 
of the N/C ratio already during the MS phase. In fact 
most of the increase in N/C is in general built during the 
MS phase, and this results from the steeper fi-gradients 
and greater compactness. The relative growths of the N/C 
ratio do not change very much from models with Z = 0.02 
to models with Z = 0.004, however there is an increase 
by two orders of a magnitude for Z = 10^^. Of course 
this large N/C enhancement is accompanied by a small 
enrichment in helium at the surface, typically of a few 



hundredths as shown by Table ||. Fig. illustrates the 
fact we find throughout this work, i.e. that the various 
effects of rotation on the internal structure, the surface 
composition and the yields are in general much higher at 
lower metallicities. 

In Fig. |l^, we notice for Z = 10^^ an increase by a 
factor 4.5 of the N/C ratio for an increase of 100 km s~^ 
of the initial rotation. As illustrated by Table H, during 
the He-burning phase the rotation velocities become all 
the same whatever the initial rotational velocities. Thus, 
in the He-burning phase we may have very different sur- 
face chemical compositions for actually similar rotation 
velocities. This is likely true for all stellar masses where 
fast rotation is present, but the effect is in general larger 
for larger masses. 

Curiously enough, at very low Z the fast rotating stars 
of intermediate masses which reach the TP-AGB phase 
(this occurs for M < 7M©) get a higher Z during this 
phase due to their enrichment in CNO elements. As an 
example, a 7 M© has a X(CNO) = 3.1 -IQ-^ which is 430 
times the initial CNO content. Thus very low Z stars may 
become higher Z stars near the end of their evolution. This 
might also affect the composition of planetary nebulae in 
low Z regions. 

8. Primary nitrogen production 

8.1. Brief recalls on the nitrogen synthesis 

Nitrogen is mainly produced in the CN branch of the 
CNO cycles within H-burning stellar zones (see Clayton 
1983; Arnett |1996 ). More precisely, three reactions occur 
to transform i^^Cto "N: ^^C{p, 7)13N(/3+, j/)"C(p, j)^'^'N. 
Nitrogen can also be produced in the ON cycle by trans- 
formation of -'^^O, but at a much slower rate. The reaction 
^^N(p, 7)^^0 which depletes nitrogen has a relatively low 
cross section enabling ^''N to accumulate with time. Thus, 
^^N is usually the daughter element, hence a secondary el- 
ement, of the CNO initially present in stars. 



Nitrogen is said to be primary (Talbot & Arnett 1974), 
if it is formed in a star not from the initial CNO ele- 
ments, but from the hydrogen and helium. Of course, the 
reactions forming primary nitrogen are those mentioned 
above, but the sequence of events is different: the forma- 
tion of primary nitrogen implies firstly the synthesis of 
some ^^C by the 3a-reaction in a helium burning region, 
then this new '^^C needs to be transported in an hydrogen 
burning region, where the CNO cycle will convert it to 
nitrogen. Thus, primary nitrogen is likely to be formed in 
stars with a He-burning core and a CNO burning shell, 
provided there is some transport mechanism between the 
two. The absence of such transport is the main reason why 
current models do not produce in general any primary ni- 
trogen. 

If ^^N is of primary origin, the -'^^N-abundance is pro- 
portional to that of the other primary heavy elements. 
While if nitrogen is secondary, the increase in the abun- 
dance of ^^N is proportional to the initial CNO content 
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and thus in the chemical history of a galaxy the ^^N- 
content will be proportional to the square of the CNO 
and metal content. These different behaviors provide the 
basic test for ascertain the primary or secondary origin of 
nitrogen. 

The observations point toward the need of primary ni- 
trogen sources at low metallicities (see below). The main 
problem is that the stellar models, unless ad hoc hypoth- 
esis are made, do not currently predict any primary ni- 
trogen. This suggests that some physical process may be 
missing in the stellar models. 

There are two other related problems. At solar metal- 
licities, the observations do not suggest the production 
of primary nitrogen. Thus, a global question is how is 
changing the respective efficiencies of the primary and sec- 
ondary productions during the evolution of galaxies. 
The other question concerns the relative importance of 
massive and intermediate mass stars in the production of 
primary and secondary nitrogen. This point is important 
in relation for the interpretation in terms of the star for- 
mation history the N/0 ratios observed at high redshifts 
(Pettini et al. |1995| ; Lu et al. |1998D . Nitro gen, primary and 
secondary, is produced in the longest and main evolution- 
ary phases. As long as its production is not well under- 
stood, we may doubt of the correctness of the models for 
these main phases of stellar evolution. 



8.2. The observations of the N/0 ratio 



There are several kinds of evidences in favour of pri- 
mary nitrogen in the early phases of the evolution of galax- 
ies. 

-1. An indication of primary nitrogen is provided by 
the study of the N/0 ratio in low metallicity stars of 
the galactic halo. The "discovery" of primary nitrogen 
was made by Edmunds & Pagel ( 1978 ) in a study of the 
N/0 ratio in such stars and in some external galaxies. 
Following this work, others authors (Barbuy 198S ; Tomkin 
& Lambert Il984l Matteucci |l986|; Carbon et al. |l987; 



see whether this sum is higher than the initial local CNO 
content of these galaxies. 

-2. A very compelling evidence for primary ^^N is pro- 
vided by the study of the N/0 ratios in ionized HII regions 
of blue compact dwarf galaxies (Thuan & Izotov 1995; 



Matteucci 1986; Carbon et al 
Henry et al. 200C| ) have shown that the ratio N/0 of nitro- 
gen to oxygen remains constant with a plateau at log N/0 
~ -1.7 in the early evolution of the Galaxy, thus imply- 
ing a primary origin of nitrogen. The limit in metallicity 
above which secondary production of nitrogen becomes 
important is difficult to fix with precision, since the tran- 
sition is progressive. It is around 12-1- log 0/H = 7.8 to 8.2 
according to Henry et al. (200C; cf. also Izotov & Thuan 
1999| ). Since for the Sun one has 12 + log 0/H =8.9, this 
means at a metallicity Z equal to Zq/12 to Zq/5. Above 
this limit, the N/C and N/0 ratios grow rapidly, implying 
that nitrogen is essentially a secondary element. It is not 
known whether the primary production stops completely 
for Z values higher than the above limit. A good way to 
check it would be to measure the sum of CNO elements 
in planetary nebulae of the SMC, LMC and Galaxy, to 



Kobulnicky & Skillman 199e; Izotov & Thuan 1999; Izotov 



& Thuan 2000). These HII regions also show a plateau of 
N/0 at log N/0 ~ -1.7 below 12 -t- log 0/H ~ 8.0, while 
above this limit the N/0 ratio is a steeply growing func- 
tion of 0/H, as due to the secondary production of ni- 
trogen. An example of a low metallicity galaxy is IZw 18, 
which has the lowest known metallicity (1/50 of solar), and 
which shows indicatio ns of p rimary nitrogen (Kunth et al. 
1995| ; Izotov & Thuan |l999|) . The stu dy of the N/0 ratio 



in spiral galaxies by van Zee et al. ( |1998[ ) well confirms 
the same result, with the difference that the authors find 
a plateau below 12 + log 0/H = 8.45, i.e. for abundance 
of heavy elements less than 1/3 solar. 

A problem was that some low Z damped Lya sys- 
tems have N/0 ratios lower than those observed in the 
HII regions of blue compact dwarf galaxies of the same 
Z (Pettini et al. |1995 ). However, the apparent discrep- 
ancy has been resolved by models of damped Lya systems 
which account for both ionized and neutral regions (Izotov 
et al. |2001| ). 

-3. Another argument for primary nitrogen production 
comes from the observed gradient of N/0 in spiral galax- 
ies. If nitrogen is purely a secondary element, the N/0 
gradient should be identical to that of 0/H. In general, 
the N/0 gradients tend to be shallower than the 0/H gra- 
dients (Vilchez & Esteban 1996). The various data on the 
galac tic gr adients of N/0 (Rud olph et al. 1997| ; G arnett 
et al. |1997| ; Ferguson et al. |l99^ ; Henry & Worthey |l999| ) 
generally show that the N/0 gradient is relatively flat at 
low metallicity Z , which supports the conclusion that the 
production of nitrogen is dominated by primary processes 
at low Z , while at solar or higher Z the similarity of the 
N/0 and 0/H gradients support the view that nitrogen 
is secondary. 

The situation is rather confused concerning the masses 
of the stars responsible for the injection of primary nitro- 
gen. There are authors supporting the o rigin of primary 
nitro gen in massive sta rs (M atteucci 1986 ; Th uan fc Izotov 
|1995| ; Izotov & Thuan |1999| ; Izotov & Thuan |2000D . Their 
main argument is the low scatter of the observed N/0 ra- 
tios at low Z. Indeed, if nitrogen is synthesized in massive 
stars, there is no time delay between the injection of ni- 
trogen and oxygen and thus a rather small scatter would 
result. On the contrary, if the primary nitrogen is made 
in intermediate mass stars, the N/0 ratio increases with 
time, since these stars release their nitrogen much later 
than massive stars do eject their oxygen. This would lead 
to a larger scatter in the observations, because galaxies 
are observ ed at various stages of their evolution. Izotov 
& Thuan (1999) suggest also that because of the inter- 
mediate mass star delay, the faster evolving massive stars 
must be a significant source of primary nitrogen in or- 
der to raise the log (N/0) ratio to the observed plateau 
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at 12+lg(0/H) ^ 7.2, a mctallicity they assume to corre- 
spond to a galactic age too short to allow nitrogen ejection 
by the intermediate masses. 

The situation may be not so clear, because some stud- 
ies found that a significant scatter does exist (Garnett 
1990| ; Skillman et al. |1997D . Also, Henry et al. ( |2000| ) have 



There is an extensive literature on AGB star models 



(see for example Forestini & Charbonnel 1997; Boothroyd 
& Sackmann 



1999 



calculated chemical evolution models which support the 
view that intermediate mass stars between 4 and 8 Mq, 
with an age of about 250 Myr, are likely to dominate the 
nitrogen production. 

What can we deduce if we accept the fact that the 
N/0 values show a great scatter at fixed value of 0/H ? 
A possibility might be that the observed scatter occurs 
because we are observing a large sample of HII regions in 
various stages of oxygen and nitrogen enrichments. This 
picture implies that most observed data should have rel- 
atively high N/0 values with fewer points, representing 
those objects experiencing sudden oxygen enrichment, lo- 
cated below the bulk of data, since presumably bursts are 
followed by relatively long periods of quiescence, with rel- 
atively higher N/0 ratios. 

The reality looks different. The distribution of points 
in N/0 vs. 0/H plane reveals that most points seem to 
be clustered at relatively low values. This suggests that 
the "equilibrium" or unperturbed locus where most HII 
regions reside is the low N/0 envelope. Thus, this sug- 
gests that the excursions caused by sudden injections of 
material are actually upward, toward the region of fewer 
points. This picture seems consistent with the lack of ev- 
idence for localized oxygen contamination fr om m assive 
stars in H II regions (Kobulnicky & Skillman 1997). The 
falloff in points above the N/0 envelope is more consis- 
tent with injections of nitrogen rather than oxygen. In 
this case, the nitrogen source might be WR stars or lu- 
minous blue variable st ars, b oth of which were considered 
by Kobulnicky et al. (1997) in their study of nitrogen- 



enriched H II regions in NGC 5253. They expect also a 
simultaneous enrichment in helium, and thus H II regions 
exhibiting high values of N/0 should also be checked for 
evidence of helium enrichment. 



8.3. The existing stellar models 

For massive stars, there is at present no model produc- 
ing primary nitrogen unless some ad hoc assumptions are 
made in order to reproduce the observed N/0 ratio at 
low metallicity (Timmes et al. 1995). In these ad hoc 
models, some mixing is permitted between the helium- 
and hydrogen-burning zones. Some primary nitrogen may 
also be produced in low-metallicity massive stars via some 
adjusted convective overshoot (Woosley & Weaver 1995] ). 
Without any physical explanation, it is difficult to under- 
stand why the production of primary nitrogen only occurs 
at low metallicities. Models of metal free Population HI 
stars (Umeda et al. 200C) produce some primary nitro- 
gen, but in too low quantitie s to re produce the observed 
plateau (see also Heger et al. 2000b ). 



Marigo |1998| , |2001| ). Up to the phase 
of thermal pulses on the AGB branch, the intermediate 
mass star models predict no primary nitrogen production. 
Only when the star enters the phase of thermal pulses, 
some He-burning products may be transported into the 
H-burning shell, thus producing some primary nitrogen. 
These models are complex and require a lot of comput- 
ing time. This is why the AGB models (Renzini & Voli 



1981; Marigo 199S) are "synthetic" models, which means 



that the model parameters follow some analytical rela- 
tions that have generally been fitted to the observations 
(this is the case for example for the minimum stellar mass 
experiencing the third dredge-up). In addition the depen- 
dence of this minimum mass on metallicity is based on 
observation. The same kind of adjustments are made for 
the occurrence of the hot bottom burning. While this may 
be useful for some purposes, it cannot be claimed that it 
represents consistent physics leading to primary nitrogen 
production. The published stellar yields for intermediate 
mass stars are generally based on such synthetic models. 
According to the models of Marigo (2001), the primary 
nitrogen production depends heavily on the parameters 
describing the hot bottom burning and the third dredge- 
up, both processes which are not adequately described in 
complete stellar models. This means that for intermediate 
mass stars the primary nitrogen production is not a fully 
consistent output. 

For completness, we also mention here that some ex- 
planations of the N/0 ratios advocates galactic processes, 
such as differential outflows of the chemical elements pro- 
duced by galactic winds (Edmunds |1990 ). Oxygen is pre- 
dominantly made in high-mass stars that undergo more 
violent explosions than intermediate mass stars, thus is 
more likely to be removed from the galaxy. This differ- 
ential outflow results in a decrease in the effective yield 
for oxygen with a corresponding increase in the N/0 
ratio. By observing massive spiral galaxies, van Zee et 
al. ( 1998| ) tried to minimize the complicating effects of 
gas outflow/inflow. They performed nitrogen and oxygen 
abundance measurements for 185 H II regions spanning a 
range of radii in 13 spiral galaxies and obtained for the 
N/0 ratios the same behavior as in low-metallicity dwarf 
galaxies. This result suggests that the observed trend in 
dwarf galaxies is not due to the outflow of enriched ma- 
terial in a shallow gravitational potential. They conclude 
that low-metallicity H II regions in all types of galaxies 
do show evidence of primary nitrogen production. 

8.4. The physics of the production of primary nitrogen 
in low Z rotating models 

We need to look with some details the physics which de- 
termines the synthesis of primary nitrogen and more gen- 
erally the particular yields at low Z. Some e ffects have 
already been examined by Meynet & Maeder ( 2001 ). We 
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Fig. 13. Variations as a function of the langrangian mass coordinate Mr of the abundances of various elements inside 
a non-rotating 3 M0 model at the metallicity Z — 10^^. Panel a) shows the chemical structure at the end of the core 
H-burning phase. Panels b) and c) at the middle and at the end of the core He-burning phase. The structure after 
one pulse along the Thermal Pulse-AGB phase is shown on panel d). 



organize the discussion in a systematic way: 

- 1. Effects of rotation in a 3 Mq model at Z = 10~^. 

- 2. Same problem at Z = 0.004 and 0.02. 

- 3. Effects of rotation in a 20 Mq at Z = 10"^ 0.004 

and 0.020. 

- 1. Figs. [1^ and |lj compare the variations of the 
abundances inside a non-rotating and a rotating 3 Mq 
model with Z — 10"^ at various evolutionary stages. At 
the end of the H-burning phase (panels a) , we notice the 



milder /x-gradient at the very edge of the core, this con- 
tributes to make slightly larger He-cores in rotating mod- 
els. This characteristic is generally larger in larger masses, 
it also persists and increases in further stages. We notice 
a significant diffusion of He and N throughout the star 
in the rotating models. At the middle and at the end of 
the He-burning phase (panels b and c), the differences in 
the chemical profiles are striking. In the non-rotating case, 
there is no new ^^C outside the convective core, and there- 
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Fig. 14. Same as Fig. for a rotating 3 Mq model at the metallicity Z = 10^^. The initial velocity on the ZAMS is 
300 km s~^, which corresponds to an average surface equatorial velocity during the Main Sequence equal to ~230 km 
s~^. Panel a) shows the chemical structure at the end of the core H-burning phase. Panels b) and c) at the middle 
and at the end of the core He-burning phase. The structure after the first five pulses along the Thermal Pulsc-AGB 
phase is shown on panel d). 



fore there is no primary -'^ N produced. While in the rotat- 
ing model, "'^^C (together with some ^^O ) is diffusing out 
the He-burning core and when it reaches the H-burning 
shell, it is turned by the CNO-cycle into primary -'^'^N, 
producing a big bump of ^"^N and a smaller one in ^^C. 
As the H-shell migrates toward the exterior, the bumps 
of primary -'^^N and ^^C also progressively exteirds toward 
the exterior. The height of these two bumps is growing 



during the He-burning phase, since diffusion is bringing 
more and more ^^C which is turned to ^^N, this explains 
the growth of the peak in ^''N at the outer edge of the 
intershell region. The abundance of ^''N in the rest of the 
intershell region is also growing with time and this is likely 
due to the inward diffusion of nitrogen from the peak. 
Typically, at the end of the He-burning phase, the ^''N- 
abundance in the intershell zone has increased by 3 orders 
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Fig. 15. Variation of the abundances of some elements in the intershell region of 3 Mq models at the end of the 
He-burning phase for the metallicities Z = 0.004 and 0.020. The initial velocities Uini are indicated. 



of a magnitude with respect to ^ N in the corresponding 
non-rotating model. At this stage, the integrated quantity 
of new nitrogen synthesized is 3.22 ■ 10"'^ Mq, while it is 
only 3.32 • 10"^ Mq in the non-rotating model. 

In further stages, the He-burning shell progresses out- 
wards, letting a degenerate CO core behind it and trans- 
forming the ^''N into ^^Ne (panels d in Figs. ^ and 
Also in the TP-AGB phase, the content in nitrogen of 
the outer convective zone increases a lot, due to both the 
facts that the outer convection zone deepens in mass and 
that the diffusion at the base of the convective zone con- 



tinues to proceed. This is saving from destruction a large 
fraction of the primary nitrogen produced earlier. At the 
stage shown in Fig. |l^ nearly the whole quantity of ^"^N is 
in the outer convective envelope. The integrated quantity 
of primary ^^N at this stage is 1.58 10"'^ Mq, i.e. about 
50% of what was present in panel c) at the end of the He- 
burning phase. This fraction of 50% does not change very 
much during the end of the TP-AGB phase, because the 
edge of the CO-core and the outer envelope stay very close 
in lagrangian coordinates. This nitrogen will be ejected by 
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the AGB star either by the superwinds or in the planetary 
nebula. 

In the corresponding non-rotating model (see Fig. |l^, 
panel d), we notice a very similar final structure with a 
large CO core surrounded by a convective envelope and 
two thin shells at the basis of it. As only differences, we 
notice the much smaller ^^N and ^^O abundances in the 
envelope. Also, because central degeneracy is higher and 
thus there is more cooling by neutrinos, the nuclear re- 
action ^^C(a,7)^^0 proceeds farther in the outer core re- 
gions than in the inner regions, leading to a kick in ^^C 
and a bump in ^^O as observed in panel d) of Fig. [l3| . 

The production of ^^C in the outer half of the zone 
between the He-core and the H-burning shell during the 
He-burning phase is relevant for the nucleosynthesis of "s- 
process" elements, since ^'^C is an efficient neutron source. 
The "s-elements" are produced when ^'^C is reached by 
the outer progression of the He-burning shell. Clearly, the 
formation of "s-elements" is strongly favoured in rotating 



stars (see also Langer et al. 1999). Amazingly, in the con- 
vective envelope, the mass fraction of the CNO elements 
is about 100 times the initial mass fraction of the heavy 
elements ! 

-2. It is interesting to compare the above results with 
those of models of a 3 Mq at higher metallicities. Fig. [T^ 
shows models with and without rotation for Z — 0.004 
and 0.020 at the end of the helium burning phase, i.e. 
corresponding to panel c) in Figs. |l^ and As usual, 
the models with zero rotation show flat curves separated 
by steep transitions due to intermediate convective zones. 
The models with rotation a.t Z = 0.004 and 0.020 show 
both the typical internal diffusion profile for ^^C and ^^O 
outside the core and noticeably the distributions of ^^N 
are rather similar to that observed for the models at Z 
= 10~^. There is still some quantity of primary nitrogen 
in the intershell model of the Z = 0.004 model, but it is 
relatively negligible in the Z = 0.02 model. The maximum 
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values of ^''N in the interior are similar in the three 3 M0 
models considered, independently of Z . 

What are the reasons of this relative constancy ? For 
the models illustrated, the central T are about the same, 
as normal for He, C, O cores of about the same mass at 
the middle of the He-burning phase in rotating models. 
However, the temperatures at the basis of the H~burning 
shells are different; logT = 7.573 at Z = 10"^ and 7.446 
at Z — 0.02. This is consistent with the fact that the Z 
= 10~^ model is much brighter (log L/Lq = 2.787 com- 
pared to 2.013 at Z = 0.02), because of its much lower 
opacity. The nuclear energy production (mainly of the 
H-shell) necessary to supply the stellar luminosity is ad- 
justed, as usual, by the temperature of the shell and not 
by the content in ^''N. The similarity of the distributions 
of "'^^N at low metallicity essentially results from the rota- 
tional transport of material from the core. We have seen 
in Sect. 3 that the fi-distribution during the MS evolu- 
tionary phase is different for different Z . However, in later 
phases the /z-contrast between the dense core and the sur- 
rounding layers is about the same and this contrast deter- 
mines the ri-gradicnt and in turn the importance of the 
diffusion. Thus, the diffusion of ^^C outside the core is not 
very different in models of same mass and rotation, and 
as a consequence the same is true for -'^'^N. 

At Z = 10"^ the gradient of ^''N between the H-shell 
and the envelope is much larger than at Z = 0.02, be- 
cause the difference between the peak of ^^N in the inter- 
shell region and the cosmic abundance in the envelope is 
also much larger. This has two consequences: firstly, the 
inward progression of the outer convective zone will bring 
relatively much more ^"^N in the envelope; secondly there 
is also more diffusion of ^^N in the envelope at very low 
Z. 



Finally, we also emphasize (cf. Meynet & Maeder gOOlj) 
that in the TP-AGB phase the distance between the He- 
and the H-burning shells is much smaller in lower Z mod- 
els. This effect will certainly influence considerably the 
occurrence and properties of the relaxation oscillations. 
Also, this effect makes the transport of ^^C from the He- 
burning shell to the H-burning shell much shorter, since 
the timescale for diffusion varies with the square of the 
distance. In this respect, the smaller intershell region in 
lower Z models also favours the increase of the abundance 
of primary nitrogen in the envelope. 

In summary, the higher production of primary 
in very low Z models results mainly from the relatively 
stronger peak of primary ^''N built by rotational mixing 
in the intershell region during the He-burning phase, a 
part of which is entering the envelope during its inward 
migration and another part is brought to the envelope by 
the diffusion, which is favoured by the smaller intershell 
region during the TP-AGB phase. 

3) Let us examine the 20 M0 models at Z = 10^^, 
0.004 and 0.020. In Fig. |l| for Z = 10-^ we see for the 
rotating model at the middle of the He-burning phase 
(panel e) the same kind of diffusion profile of ^^C and ^^O 
outside the core, as in the corresponding 3 Mq model. 



A similar, although slightly smaller peak of primary ^^N 
is built between the core and the H-shell. Contrarily to 
smaller masses, where there is no central C-burning, the 
intershell region remains large. The mainly primary ^^N 
in this region will of course contribute to the yield, as well 
as that in the outer envelope. The abundance of ^'*N is 
increasing in the envelope during He-burning and later 
phases. Since here, contrarily to the low mass models, 
there is no inward migration of the envelope, the increase 
of ^"^N in the envelope is only due to the diffusion from the 
^■^N gradient in the H-burning shell. As in smaller masses, 
most of the ^"^N in the envelope is primary. 

When we do a similar study in the 20 Mq models at Z 
= 0.004, we notice that there is only a negligible amount 
of primary ^"^N produced and the abundance of ^"^N in the 
envelope of the final models is within a few percents the 
same in the rotating and non-rotating cases. 

In conclusion, we see that there is still some primary 
nitrogen produced ai Z — 0.004 in intermediate mass mod- 
els, as shown by the 3 Mq case, but nothing in the high 
mass stars. 



9. The stellar yields 

9.1. Ma, Mco 3nd the mass of the remnants, Afrcm 

The computation of the stellar yields, i.e. the quantities 
of an element newly produced by a star, necessitates, as 
a first step, the estimate of the masses Ma of the helium 
cores, Mco of the carbon-oxygen cores, and of the masses 
Mrem of the rcmnauts. These quantities are used for two 
purposes: 1) To obtain the oxygen yields using the rela- 
tion between Ma and the oxygen yields by Arnett ( 1991 ) . 
Our models only give an upper limit of the oxygen yields, 
since they are stopped at the end of the carbon or helium 
burning phase, i.e. at phases where oxygen has not yet 
been depleted in the inner regions. 2) To obtain the mass 
of the remnants using a relation between Mqo Si,nd Mrem- 
The knowledge of this quantity allows us to calculate the 
mass of the elements expelled by the supernova. In this 
work we use the same Mj-em vs Mco relation as Maeder 
(1992). For the intermediate mass stars, we have taken as 



remnant masses, the mass Mco of the CO core. 

Both relations, Mq versus the oxygen yields and 
Mco versus Mrcm are deduced from non-rotating mod- 
els. Ideally one should have used relations obtained from 
rotating models and more precisely from rotating models 
treating the effects of rotation as we did. If such compu- 
tations would have been available, would the stellar yields 
be the same ? For what concern helium, nitrogen and car- 
bon (although to less extent) the answer is yes. Indeed the 
parts of the stars which contribute the most to the yields 
in helium, nitrogen and carbon are in too far out portions 
of the star for being affected by the Mco versus Mj-em 
relation. For the oxygen and heavy elements yields the 
situation is more complicated. For the moment, we can 
say that, unless strong and very rapid mixing episodes 
take place after the end of the carbon burning phase, it 
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Fig. 17. Variation as a function of the initial mass of the 
stellar yields in ^''N for different metallicities and rota- 
tional velocities. The continuous lines refer to the models 
at Z = 10^^ of the present paper, the dotted lines show 
the yield s from the models at Z — 0.004 from Maeder & 
Meynet ( 2001 ), the dashed lines present the yields for two 
solar metallicity models (present work). The filled squares 
and circles indicate the cases without and with rotation 
respectively. In this last case Uini = 300 km s~^. The 
crosses are for the models of Woosley and Weaver (1995) 
at Z = O.IZ0, the empty squares for the yields from van 
den Hoek & Groenewegen ( 1997] ) at Z=0.004, the empty 
triangles are for solar metallicity models of van den Hoek 
& Groenewegen (1997) up to 8 M© and of Woosley and 
Weaver (1995) above. 



is likely that the yields in oxygen and heavy elements ob- 
tained here are good estimates of the yields one would 
have obtained from rotating presupernovae models. 

In Table |[ Ma, Mco and Afrcm, are indicated for dif- 
ferent initial mass stars with various metallicities and ini- 
tial rotation velocities. For initial masses inferior or equal 
to 7 Mq , the core masses are estimated after the end of the 
He-burning phase, during the TP-AGB phase. For higher 
initial mass models, Ma and Mco are evaluated at the 
end of the C-burning phase at Z = 10^^, except for the 9 
M0 models with Vmi = 0, 200 and 400 km s~^ for which 
the core masses are estimated at the beginning of the C- 
burning phase. For the models at Z = 0.004 (Maeder & 
Meynet 2001), the core masses are estimated at the end of 
the He-burning phase. We define Ma as the mass interior 
to the shell where the mass fraction of helium becomes 
superior to 0.75. For the rotating 40 and 60 M© models 
at Z = 10~^, the diffusion of helium outside the He-core 
is so great that the above definition yields an unrealis- 
tic high value for the helium core. For these models, we 




200 

mean velocity [km s~'] 

Fig. 18. Variation as a function of the mean velocity dur- 
ing the Main Sequence, v, of the stellar yields in ^^N for 
a 9 and a 20 M© model at Z = 10"^. 



choose to fix the outer border of the He-core at the po- 
sition where the hydrogen mass fraction goes to zero. Let 
us note that for the other initial masses this alternative 
definition of the He-core does not change the results pre- 
sented in Table ^. The mass of the carbon-oxygen core 
Mco is the mass interior to the shell where the sum of 
the mass fractions of ^^C and ^^O is superior to 0.75. 

From Table |^, we note, that rotation in general in- 
creases the masses of the helium and CO cores for the 
massive stars. The reason is that for higher rotation, the 
intermediate convective zone, associated to the H-burning 
shell, disappears more quickly. Since the H-burning shell 
is not replenished in hydrogen, it can more quickly mi- 
grate outwards and thus produces the general increase of 
the He-core mass. However, for the models at Z ^ 10^^ 
we notice a saturation effect in the increase of M^ for 
higher rotation and even a decrease for a very high ro- 
tation (see the 9 and 20 M© models in Table |3|). This 
behaviour results from the following opposite effect: when 
rotation increases, the diffusion becomes so efficient that 
large amounts of hydrogen are brought from the radiative 
envelope into the H-burning shell. This slows down its out- 
ward progression and thus does not produce the growth of 
the He-core mass in the He-burning phase. As to Mco, 
we see that the variations of Mco with Vini follows those 
of Ma- 
in conclusion, we find that fast rotation in general in- 
creases Mco and thus will also increase the yields in a- 
elements. This is true at very low Z, but not necessarily at 
solar metallicity, because there fast rotating massive stars 
will experience high mass loss. 
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Table 3. Masses of the helium cores, of the carbon-oxygen 
cores and of the remnants for different initial mass star 
models with and without rotation a.t Z = 10^^ and 0.004. 
The masses are in solar mass and the velocities in km s~^. 
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Table 4. Stellar yields for helium, carbon, nitrogen, oxy- 
gen and the heavy elements from different initial mass 
stars at Z = 10^"" with and without rotation. The initial 
stellar masses are in solar units. The stellar yields are in 
solar mass and the velocities in km s~^. The contributions 
of the stellar winds have been accounted for and are indi- 
cated in parenthesis when they exceed one percent of the 
total yield. 



z = 10-^ 



M 


^^ini 


mpHc4 


"ipci2 




mpoie 




60 





l.lc-t-1 


8 


4c- 1 


1.6c-4 


l.lc+1 


1.3o+l 




300 


1.0c-(-l 


6 


lc-1 


7.0C-4 


1.60+1 


1.90+1 






(0.2O-I-1) 






(0.4c-4) 






40 





6.9c+0 


6 


3c- 1 


l.Oc-4 


5.O0+O 


7 3e+0 




300 


S.Oc+O 


5 


8c- 1 


7.6C-4 


5.80+0 


8.20+0 


20 





3.1C+0 


2 


2c- 1 


4.6C-5 


5.60-1 


1.70+0 




200 


2.5C+0 


2 


8c- 1 


4.2C-5 


1.20+0 


2.80+0 




300 


3.3C+0 


2 


9c- 1 


3.6C-4 


9.90-1 


2.30+0 




400 


3.5e+0 


2 


9c- 1 


7.7C-3 


6.20-1 


2.O0+O 


15 





2.0C+0 


1 


2c- 1 


3.5C-5 


1.90-1 


8.60-1 




300 


1.9C+0 


1 


9c- 1 


4.1C-4 


3.90-1 


1.4o+0 


9 





l.Oc+0 


1 


lc-2 


2.2C-5 


4.O0-3 


5.I0-2 




200 


l.lc+0 


3 


2c-2 


7.7C-4 


4.O0-3 


2.O0-I 




300 


l.lc+0 


3 


6c- 2 


3.5C-3 


4.O0-3 


1.60-1 




400 


l.lc+0 


4 


8c- 2 


9.2C-3 


4.O0-3 


2.O0-I 


7 





6.9C-1 


1 


3c- 2 


1.6C-5 


1.20-3 


1.4o-2 




300 


8.O0-I 


1 


8c- 2 


4.2C-3 


5.60-3 


2.80-2 


5 





4.I0-I 


6 


9c- 3 


l.Oc-5 


5.60-4 


7.70-3 




300 


4.5C-1 


4 


5c-3 


4.3C-3 


1.80-3 


I.I0-2 


3 





9.O0-2 


4 


7c-4 


3.3C-6 


8.O0-6 


7.O0-4 




300 


1.40-1 


2 


5c-3 


1.6C-3 


7.O0-4 


5.I0-3 


2 





9.80-2 


1 


5c-3 


1.7C-6 


1.90-5 


1.8e-3 




300 


1.2C-1 


7 


7c- 3 


6.4c-4 


I.I0-3 


l.Oe-2 



For the intermediate mass stars, the situation is more 
complicated, since in addition to the effects just mentioned 
above, the mass of the helium core also results from the 
inward progression in mass of the outer convective zone 
during the AGB phase. The deeper this inward extension, 
the smaller is Ma. 



9.2. The stellar yields 

In Tables || and ^, the stellar yields for helium, carbon, 
nitrogen, oxygen and for the heavy elements are given for 
different initial stellar masses with various initial metallic- 
ities and rotational velocities. Except for oxygen in mas- 
sive stars, our models have reached a sufficiently advanced 
evolutionary stage for the above yields to be directly de- 
duced from our models. The quantity of an element x, 
newly produced by a star of initial mass to, i.e. the stellar 
yields in x, is given by 



A/tin 



XJdlTlr 



A/rcn 



where Mfin is the mass of the star at the end of its evo- 
lution, Xx{mr) is the mass fraction of element x at the 
langrangian mass coordinate rur inside the star and X° 
is the initial abundance of element x. We consider that 
the SN ejecta in carbon consist of the C-distribution in 
the star as it is at the end of the C-burning phase, but 
we count only the layers which are external to the mass, 
where carbon is not depleted by the further nuclear burn- 
ing stages (see Maeder 1992 ). 

We also accounted for the effec ts of the stellar winds 
on the yields as is done in Maeder ( 1992 ) . In Tables ^ and 
H, we have indicated in parenthesis the contribution of 
the winds when it exceeds one percent of the total stellar 
yields. For the metallicities considered here, the effects of 
the winds are small and only modify the yields in helium 
and nitrogen. 

In Fig. the stellar yields in ^^N are plotted as 
a function of the initial mass for different metallicities 
and rotation velocities. Our non-rotating models (filled 
squares along the continuous lines) show very small yields 
dX Z = 10~^, as expected from a pure secondary origin of 
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Fig. 19. Variation of the stellar yields in oxygen and carbon as a function of the initial mass, at three different 
metallicities. The black triangles are for the present non-rotating models, the black circles are for the present rotating 
models. The empty squares and pentagons are for the stellar yields of Maeder ( |l992| ) at the metallicity Z = 0.001 and 
0.020 respectively. 



^"^N. When the effects of rotation are included (filled circles 
along the continuous lines), the yields for the intermediate 
mass star models at Z = 10~^ become of the same order 
of magnitude as the yields for the corresponding models 
at Z = 0.004. In this mass range, the yields present thus 
a very weak metallicity dependence. At solar metallicity, 
the yields given by the rotating and non-rotating mod- 
els are identical, this well illustrates the smaller effects of 
rotation on the nitrogen yields for the higher metallicities. 

Figure |8| shows that the stellar yields at Z = 10"^ 
strongly depend on the rotational velocity. Starting from 
the point at u = 230 km s~^ for the 9 M0, an increase 
of 20% of the mean velocity on the MS, already suffices 
to double the yield in nitrogen. One notes also the strong 
increase obtained for the 20 Mq when v passes from 240 
to 325 km s^^. This means that if massive stars have sufh- 
ciently high initial velocities at low Z, they might also play 
a role in the production of primary nitrogen (see below) . 

Fig. |l^ shows the yields in carbon and oxygen for 
Z = 10-^ 0.004 and 0.020 for models with and without 



rotation. We notice that the yields of carbon and oxygen 
(and other heavy elements) are much less affected by ro- 
tation than the yields in nitrogen. In general, the yields 
in carbon, oxygen (and heavy elements) are increased by 
rotation (cf Heger 1998), this is a result of the generally 
bigger value of Mco when rotation is included, as shown 
above. At solar metallicity, rotation decreases the yield in 
oxygen because of the higher mass loss, which also pro- 
duces an increase in the amount of carbon ejected. 

How the present stellar yields compare with the yields 
from other authors ? The situation for nitrogen can be seen 
on Fig. |r^, where yields of van den Hoek & Groenewegen 
( 1997t see their tables 9 and 17) and of Woosley and 
Weaver ( [1995 , see their models S and P) are plotted. The 
yields at Z = 0.004 of van den Hoek & Groenewegen 
( |1997|) for the masses between 5 and 8 Mq are higher by 
about an order of magnitude than the yields of their lower 
initial mass models. They are nearly at the same level as 
the yields obtained from the solar metallicity models. This 
huge enhancement of their ^"^N stellar yields in this mass 
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Table 6. Integrated yields in helium, carbon, nitrogen, oxygen and heavy elements (see text). 



Z 12+log(g) ma TOu t>i„i Phs Pc Pn Po Pz Log |S- Log ^ 



10- = 


5.74 


20 


120 





4.9E-3 


3.8E-4 


7.4E-8 


4.1E-3 


5.1E-3 


-1.034 


-4.743 






8 


120 





7.8E-3 


4.9E-4 


1.3E-7 


4.2E-3 


5.8E-3 


-0.936 


-4.515 






2 


120 





l.lE-2 


5.3E-4 


2.1E-7 


4.2E-3 


5.9E-3 


-0.900 


-4.303 






2 


60 





9.2E-3 


4.2E-4 


1.8E-7 


1.9E-3 


3.7E-3 


-0.659 


-4.013 


0.004° 


8.35 


20 


120 





4.3E-3 


2.5E-4 


2.9E-5 


4.7E-3 


5.6E-3 


-1.269 


-2.205 






8 


120 





7.3E-3 


3.3E-4 


4.4E-5 


4.9E-3 


6.4E-3 


-1.175 


-2.050 






2 


120 





l.lE-2 


4.6E-4 


6.5E-5 


4.9E-3 


6.6E-3 


-1.024 


-1.878 






2 


60 





9.5E-3 


3.6E-4 


5.5E-5 


2.4E-3 


4.0E-3 


-0.824 


-1.638 


0.020*' 


8.93 


2 


120 





8.9E-3 


4.4E-4 


3.6E-4 


3.3E-3 


5.9E-3 


-0.881 


-0.963 


10" = 


5.74 


20 


120 


300 


4.8E-3 


3.3E-4 


4.1E-7 


6.4E-3 


7.7E-3 


-1.286 


-4.200 






8 


120 


300 


7.7E-3 


5.3E-4 


3.8E-6 


6.8E-3 


9.0E-3 


-1.106 


-3.245 






2 


120 


300 


l.OE-2 


6.2E-4 


l.lE-4 


7.4E-3 


9.5E-3 


-1.080 


-1.833 






2 


60 


300 


l.OE-2 


5.4E-4 


3.4E-5 


2.7E-3 


5.0E-3 


-0.698 


-1.898 


0.004° 


8.35 


20 


120 


300 


3.7E-3 


3.8E-4 


2.8E-5 


7.0E-3 


7.8E-3 


-1.265 


-2.396 






8 


120 


300 


6.9E-3 


5.7E-4 


5.8E-5 


7.4E-3 


9.3E-3 


-1.112 


-2.109 






2 


120 


300 


l.OE-2 


6.2E-4 


l.lE-4 


7.4E-3 


9.5E-3 


-1.080 


-1.833 






2 


60 


300 


8.8E-3 


4.6E-4 


9.9E-5 


3.4E-3 


6.0E-3 


-0.871 


-1.536 


0.020'' 


8.93 


2 


120 


300 


l.lE-2 


1.2E-3 


3.9E-4 


2.7E-3 


6.5E-3 


-0.370 


-0.844 



Models of Macdor & Moynct SL 
Models of Meynet & Maoder EOOO 



Table 5. Same as Table | for the metalhcity Z=0.004. 



M 


^iiii 


mpHc4 


mpci2 


Z = 0.004 

mpNi4 




mpz 


60 





9.3C+0 


6.7C-1 


6.8e-2 


1.2e-f 1 


1.4e-|-l 






(2.1O+0) 




(1.5e-2) 






40 





6.3C+0 


4.0C-1 


4.3e-2 


5.3e-|-0 


7.3e-|-0 






(0.2O+0) 




(0.2e-2) 








300 


5.2C+0 


4.9C-1 


3.9e-2 


7.3e-|-0 


l.Oe+1 






(4.1O+0) 




(2.6e-2) 






25 





3.7C+0 


l.lc-1 


2.1e-2 


1.9e-f0 


3.1e-|-0 




300 


3.1O+0 


2.4o-l 


2.2e-2 


2.6e+0 


4.7e+0 






(O.lc+0) 




(0.2e-2) 






20 





3.0C+0 


8.6C-2 


1.6e-2 


8.5e-l 


1.8e-|-0 




200 


2.7C+0 


1.5C-1 


1.7e-2 


1.3e-f0 


2.7e-|-0 




300 


2.5O+0 


2.I0-I 


1.7e-2 


1.4e-|-0 


2.9e4-0 










(O.le-2) 








400 


2.5O+0 


2.80-1 


1.7e-2 


1.4e-|-0 


3.2e4-0 










(O.le-2) 






15 





2.O0+O 


3.50-2 


l.Oe-2 


2.5e-l 


8.0e-l 




300 


2.0C+0 


1.90-1 


1.4e-2 


3.9e-l 


1.4e-|-0 


12 





1.4o+0 


8.30-2 


6.8e-3 


2.7e-2 


4.9e-l 




300 


l.Sc+O 


9.30-2 


1.6e-2 


7.0e-2 


6.0e-l 


9 





1.2C+0 


1.4o-2 


5.4e-3 


4.0e-3 


1.6e-2 




300 


1.3C+0 


4.60-2 


1.6e-2 


4.0e-3 


1.8e-l 


3 





1.60-1 


8.80-3 


l.le-3 


6.5e-4 


l.le-2 




300 


1.3e-l 


7.6e-4 


2.4e-3 


-2.7o-4 


3.6e-3 



nitrogen. The values of the parameters (minimum core 
mass for the third dredge-up, third dredge-up efficiency, 
scaling law for mass loss on the AGB, core mass at which 
the hot bottom burning is assumed to operate) have been 
chosen in order to reproduce various observational con- 
straints, as e.g. the luminosity function of carbon stars 
or the abundances observed in planetary nebulae (see van 
den Hoek & Groenewegen 1997). Interestingly, our rotat- 
ing models are well in the range of these parameterized 
yields. Thus rotation, not only naturally leads to the pro- 
duction of primary nitrogen, but also predicts yields in 
primary nitrogen at a level compatible with those deduced 
from previous parameterized studies. 

One also notes that our yields at Z— 0.004 (both 
from rotating and non-rotating models) are inbetween 
the yields at Z-O.O OA and Z=0.020 of van den Hoek & 
Groene wegen ( 1997 ) and between those of Woosley and 
Weaver (|l995| ) at Z ~ 0.002 and 0.020. This indicates that 
at higher metallicity, the present yields in nitrogen seem to 
be in agreement with the yields of other authors. Similar 
conclusion are reached when comparisons are made be- 
tween our yields in carbon and oxygen with those of these 
authors. Our carbon and oxygen yields also compare well 
with those obatined by Maeder ( 1992 ) (see Fig. 19). 



range is a consequence of accounting, in a parameterized 
way, for the effects of the third dredge-up and of the hot 
bottom burning, which enable the production of primary 



9.3. Net yields and comparison with the observations 

In order to evaluate the impact of these new yields (see 
Tables | and |) on a galactic scale, we use them in a 
very simple model of galactic chemical evolution making 
use of the closed box, instantaneous recycling approxima- 
tions and supposing a constant star formation rate. We 
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Fig. 20. Simplified model for the galactic evolution of the 
C/0 ratio as a function of the 0/H ratio (in number). The 
dashed and continuous lines show the results deduced from 
the non-rotating and the rotating models respectively (see 
text). The range of the initial masses used for computing 
the integrated yields are indicated. The empty symbols 
show the results when only stars more massive than 8 
Mq are considered. The initial velocity is indicated. The 
shaded area shows the region where most of the obser- 
vations of extragalactic HII regions and stars are located 
(see e.g. Gustafsson et al 1999; Henry et al. 200C| ). 



are of course fully aware of the roughness of these approx- 
imations, but our intention here is just to estimate the 
relative effects of rotation on the chemical yields. 

In the conditions of this simplified chemical evolution 
model, the ratio Xi jxj of the mass fractions of the elements 
Xi and Xj in the interstellar medium are given by ^ = 

ffi", where yi and yj are representative time-independent 

approximations of the integrated yields of the elements i 
and J from a stellar generation. The integrated yield of an 
element a;, P^-, is defined as the mass fraction of all stars 
formed, which is eventually expelled under the form of the 
newly synthesized element x: 



Pr. 



where <l>(m) is the Initial Mass Function (IMF). Here we 
choose a Salpeter IMF, normalized so that 

"120Mq 

^{m)dm — 1, 

O.1M0 

The masses md and niu limit the mass range of the stars 
having contributed to the chemical evolution of the inter- 
stellar medium at the epoch considered. 



Fig. 21. Same as Fig. |g for the N/0 ratio. 



In Table ^, the integrated yields are given for various 
metallicities, initial rotational velocities and values of md 
and TOu. The integrated yields for the solar metallicity 
have been deduced from the models of Meynet & Maeder 
( 2000| ) . These models have been computed with a different 
shear diffusion coefficient, and for a different prescription 
of the mass loss rate than the present models. Therefore 
they do not belong to the homogeneous set of data consti- 
tuted by the models at and 0.004. Despite these 
differences in the physical ingredients, their properties are 
well in the lines of the results obtained at lower metal- 
licity. This is the reason why, we have complemented the 
data of Table || with the integrated yields obtained from 
these models at solar metallicity. 

In Figs 1^ and the C/0 and N/0 ratios, obtained 
by simply taking the ratios of the corresponding integrated 
yields, are plotted as a function of 0/H. To disentangle 
the still controversial role of the intermediate mass stars 
and the effects of rotation, we have taken several values 
of the upper and lower mass limits. We have considered 
at Z = 10"^ the 20 to 120 Mq interval, that of 8 to 120 
Mq and that of 2 to 60 Mq. In each case, the models 
without rotation and with an initial rotation of 300 km/s 
(i.e. an average of 230 km/s during the MS phase) have 
been considered. In addition, for the case of 8 to 120 Mq 
an initial rotation of 400 km/s has also been considered. At 
Z — 10^^, we notice that rotation only slightly decreases 
the C/0 ratio, the effect is a bit larger when only massive 
stars are considered. This behaviour is due to the growth 
of Mco with rotation, which favours the production of 
O more than that of C. However, we emphasize that the 
effect of the mass interval is more important than rotation. 
When the intermediate mass stars are included, the C/O 
ratio is as expected much larger. In summary, at low Z, 
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the diagram C/0 vs. 0/H is particularly sensitive to the 
mass interval. 

Between Z = 0.004 and 0.02, the main effect influenc- 
ing the C/0 ratio is no longer the value of Mco as above, 
but the effects of stellar winds and their enhancement by 
rotation. Rotation favours a large C/0 ratio, because ro- 
tating models enter at an earlier stage into the Wolf-Rayet 
phase than their non-rotating counterparts. As a conse- 
quence, in rotating models at Z — 0.02, great quantities 
of carbon are ejected by the massive stars through their 
stellar winds, when they become a Wolf-Rayet star of the 
WC typ e. Thi s is quite in agreement with the results by 
Maeder ( 1992 ) who showed that when the mass loss rates 
are high, most of the carbon is produced and ejected by 
massive stars through their stellar winds. In non-rotating 
models, the new mass loss rates used here are much lower 
than those used in 1992 because now the mass loss rates 
account for the clumping effects in the Wolf-Rayet stellar 
winds. 

The diagram N/0 vs 0/H (Fig. ^) has a different sen- 
sitivity to the mass interval and rotation. At very low Z, 
we notice a very high sensitivity to rotation when the low- 
est mass limit is at 2 or 8 M0 ; this is due to the production 
of primary nitrogen. The increase in the N/0 ratio may 
reach more than 2 orders of a magnitude. Contrarily to the 
previous diagram, the N/0 ratio is not sensitive at all to 
the mass interval for models without rotation. Thus, the 
combination of the diagrams C/0 vs 0/H, more sensitive 
to the mass interval, and of the diagram N/0 vs. 0/H, 
more sensitive to rotation, may be particularly powerful 
to disentangle the two effects of rotation and mass inter- 
val, and to precise the properties of the star populations 
responsible for the early chemical evolution of galaxies. 

At the present stage, when we compare our results to 
the observations we may derive the following tentative 
conclusions, which could change if the data further im- 
prove. The C/0 vs 0/H diagram does not seem favorable 
to enrichments by only very massive stars in the range 20 
to 120 Mq; contributions from stars down to 8 or 2 Mq 
may be needed, depending on the exact slope observed 
at low Z in Fig. As to the N/0 vs 0/H diagram, no 
model without rotation is able to account for the observed 
plateau, moreover contributions from only stars above 20 
Mq seem difficult. The observed plateau at log N/0 = 
-1.7 strongly supports rotating models including the large 
contribution from intermediate mass stars down to, either 
2 Mq if these stars have the same average rotation as 
in Pop. I stars, or down to only 8 Mq if the rotations are 
faster as suggested by Maeder et al. ( 1999] ). In this respect, 
it would be sufficient that the average rotational velocities 
during the MS phase are larger by about 80 km/s. 

We must temperate these conclusions by the follow- 
ing remark (cf. also Meynet & Maeder |2001| ), related to 
a current problem in the chemical evolution of galaxies. 
Nitrogen is ejected mainly by AGB stars with ejection ve- 
locities of a few 100 km s~^, while oxygen is ejected by su- 
pernovae at much higher velocities of 10^ km s~^ or more. 
Thus, a fraction of the oxygen produced may escape from 



the parent galaxy, leading to a higher N/0 ratio than in 
the simple estimate made here. 

10. Conclusions 

We have investigated with some details the physics of very 
low Z models with Z = 10^^ . Even if we have little chance 
before long to observe such star populations, these mod- 
els are most relevant for the early chemical evolution of 
galaxies. The models enable us to calculate the yields in 
CNO, He and heavy elements necessary to study the early 
evolution of galaxies. The diagrams showing log N/0 vs 
log 0/H and log C/0 vs 0/H are particularly powerful to 
infer properties of the early star generations in galaxies. 

This work shows the large differences brought by the 
present models with rotation. In future we shall calculate 
the detailed yields of the various heavy elements. 
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